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In this paper we consider the one-dimensional Hubbard model and study the deviations from 
the ground state values of double occupation which result from creation or annihilation of holons, 
spinons, and pseudoparticles. These quantum objects are such that all energy eigenstates are de- 
scribed by their occupancy configurations. The band-momentum dependence of the obtained double- 
occupation spectra provides important information on the degree of localization/delocalization of 
the real-space lattice electron site distribution configurations associated with the pseudoparticles. 
We also study the band-momentum, on-site electronic repulsion, and electronic density dependence 
of the pseudoparticle energy bands. The shape of these bands plays an important role in the 
finite-energy spectral properties of the model. Such a shape defines the form of the lines in the 
momentum-energy /frequency plane where the peaks and edges of the one-electron and two-electron 
spectral weight of physical operators are located. Our findings are useful for the study of the 
one-electron and two-electron spectral weight distribution of physical operators. 
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I. INTRODUCTION 

Recently, there has been a renewed experimental interest in the exotic one-electron and two-electron spectral 
properties of quasi-one-dimensional materials Some of 

these experimental studies observed unusual finite-energy/frequency spectral properties, which are far from being well 
understood. For values of the excitation energy larger than the transfer integrals associated with electron hopping 
between the chains, the one-dimensional (ID) Hubbard model [2ll I22I I23L I24I |25| is expected to provide a good 
description of the physics of these materials 0, H 0, 0] . Moreover, recent angle-resolved ultraviolet photoemission 
spectroscopy revealed very similar spectral fingerprints from both high-T c superconductors and quasi-one-dimensional 
compounds @ . The similarity of the ultraviolet data for these two different systems could be evidence of the occurrence 
of a charge-spin separation associated with holons and spinons. The anomalous temperature dependence of the spectral 
function could also indicate a dimensional crossover 

U HI El EH. The results of Refs. |2£j also suggest that 
the unconventional spectral properties observed in two-dimensional (2D) materials could have a ID origin. Thus the 
holons and spinons could play an important role in spectral properties of both ID and 2D low-dimensional materials. 

The present study is related the to both the investigations of Refs. [3lJ and [13 • In reference [3l| the non- 
perturbative organization of the N electrons and N h — [2N a —N] electronic holes, which results from the many-electron 
interactions, was studied. In that paper it was found that all energy eigenstates of the model can be described in 
terms of occupancy configurations of a c = ±1/2 holons, a s — ±1/2 spinons, and c pseudoparticles where a c and a s 
are fy-spin and spin projections respectively. Moreover, the c, v pseudoparticles (and s, v pseudoparticles) associated 
with Takahasi's charge (and spin) ideal string excitations of length v |22ll31| were found to be 77-spin singlet 2z/- holon 
(and spin singlet 2^-spinon) composite quantum objects. (The c and a, v pseudoparticles are obtained naturally and 
directly from analysis of the Bethe-ansatz solution [3l|.) In that reference, the ±1/2 holons (and ±1/2 spinons) 
which are not part of such 2^-holon composite c, v pseudoparticles (and 2^-spinon composite s, v pseudoparticles) 
were called ±1/2 Yang holons (and ±1/2 HL spinons). In the designation Yang holon (and HL spinon) Yang refers to 
C. N. Yang (HL stands for Heilmann and Lieb) j^l). Following the notation of Ref. [31j, here we call N c and N a> „ the 
numbers of c pseudoparticles and a, v pseudoparticles belonging to branches a = c, s and u — 1,2, 3, ... respectively. 
On the other hand, it was found in Ref. [3l| that the charge excitations associated with the c pseudoparticles are 
independent of the 77-spin-holon and spin-spinon excitations. Moreover, the results of Ref. [32| reveal that the above 
quantum objects can be expressed in terms of rotated electrons. According to these results the first step of the non- 
perturbative diagonalization of the model corresponds to a mere unitary rotation in Hilbert space which maps the 
electrons onto rotated electrons. For rotated electrons double occupation is a good quantum number for all values of 
the on-site repulsion. Interestingly, such a rotation corresponds to a unitary transformation previously introduced by 
Harris and Lange |33|. Such a transformation is shown in Ref. |32l | to play a key role in the expression of the holon 
and spinon number operators in terms of electronic operators for the whole parameter space of the model. The N c 
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c pseudoparticles (and N c spinons) describe the charge part (and spin part) of the N c rotated electrons which singly 
occupy lattice sites. Thus the charge and spin degrees of freedom of such a 7V c -rotated-electron fluid separate. For 
each c pseudoparticle there is a chargeon [^U and a rotated-electronic hole. Such a chargeon corresponds to the charge 
part of a rotated electron which singly occupies a site. We note that the numbers iV c , [N — N c ]/2, and [N h — N c ]/2 
equal the number of rotated-electron singly occupied sites, doubly occupied sites, and empty sites respectively. The 
non-perturbative organization of the electronic degrees of freedom also includes a [N — iV c ]-rotated-electron fluid 
associated with the rotated-electron pairs occupying [N — N c ]/2 sites. Each such a doubly occupied site corresponds 
to a spin singlet rotated electron pair which is nothing but a —1/2 holon. Such a two-fluid picture shows some formal 
analogies to those which characterize a superconductor or a superfiuid. Finally, the rotated-electron [N h — N c ]/2 
empty sites are nothing but the +1/2 holons. 

The main goal of this paper is the study of the c pseudoparticle, a, v pseudoparticle, —1/2 Yang holon, and —1/2 HL 
spinon double-occupation spectra. These functions are evaluated by use of the Hellmann-Feynman theorem [34 , l35| . 
In the case of the pseudoparticles we find band-momentum dependent double-occupation spectra. Analysis of these 
spectra provides important information about the localization/deslocalization degree of the real-space lattice electronic 
occupancy configurations which describe the pseudoparticles. We find that for some of the pseudoparticle branches 
such a degree of localization/delocalization is strongly dependent on the value of the pseudoparticle band momentum. 
Our analysis leads to a double-occupation functional whose coefficients are the double-occupation spectra. This study 
confirms that the complex electron site distribution configurations of the real-space lattice which describe the energy 
eigenstates are dependent on the value of U/t. We also consider U/t — > oo electron double-occupation selection rules. 
These rules are used in Ref. in finding exact rotated-electron selection rules which provide for all values of the 
on-site electronic repulsion the number of holons and spinons of the final states which contribute most significantly 
to few-electron correlation functions. These selection rules imply that a large number of matrix elements between 
energy eigenstates are vanishing. This simplifies the derivation of general expressions for correlation functions at finite 
excitation energy. This problem is studied in Refs. [36l l37| . 

In addition, we study the band-momentum, on-site electronic repulsion U, and electronic density n dependence of 
the pseudoparticle energy bands. Such pseudoparticle energy bands are the coefficients of the first-order term of the 
finite-energy Landau-liquid functional. That functional was obtained in Ref. |38|| by expanding the c pseudoparticle, 
c, v pseudoparticle, and s, v pseudoparticle band-momentum distribution functions around their ground-state values. 
This is similar to the case of the quasi-particle energy functional of Fermi liquid theory |39i |40j . The pseudoparticle 
bands are the basic blocks of the energy spectra of the elementary excitations of the ID Hubbard model, as for 
example the ones studied in Refs. [i H W A li a EL EH li^ . 

The finite-energy theory constructed here and in Refs. |3lll32l l36l l37| is applied elsewhere to the study of the finite- 
energy /frequency one-electron and two-electron spectral properties of quasi-one dimensional materials. Reference 01 
is an application of a preliminary version of our finite-energy theory. The evaluation of the spectral functions of the 
model reveals that the one-electron and two-electron spectral weight of physical operators contains peaks and/or edges 
which are located on well-defined lines in the momentum k and energy/frequency to plane. Importantly, the shape 
of these peak and edge lines is fully defined by the band-momentum dependence of the pseudoparticle energy bands 
studied in this paper. Moreover, the expressions of the finite-energy correlation functions involve these pseudoparticle 
bands Is?!. Therefore, our investigation is important for the study of the one-electron and two-electron spectral- 
weight distribution of physical operators. The study of the localization-delocalization degree of the real-space lattice 
electronic occupancy configurations which describe the pseudoparticles also contributes to important information 
about the physics behind these elementary quantum objects. Such a study also provides information about the 
electron site distribution configurations of the real-space lattice which describe the energy eigenstates. 

The paper is organized as follows: In Sec. II we introduce the ID Hubbard model. The concepts of CPHS ensemble 
space and subspace are introduced in Sec. III. Here CPHS stands for c pseudoparticle, holon, spinon. In Sec. IV the 
double occupation functional and its pseudoparticle, —1/2 Yang holon, and —1/2 HL spinon spectra are obtained and 
discussed. In this section we also consider U/t — > oo selection rules in the values of double occupation deviations which 
result from A/"-electron excitations. In Sec. V we study the band-momentum dependence of the energy pseudoparticle 
bands. Finally, the concluding remarks are presented in Sec. VI. 



II. THE ID HUBBARD MODEL 



In a chemical potential \i and magnetic field H the ID Hubbard Hamiltonian can be written as, 




(1) 



a—c, s 
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where the Hamiltonian, 



H 



50(4) 



H 



H 



(17/2) N + (17/4) N a . 



(2) 



has 50(4) symmetry [43, |4j| and 



H H =T + UD, 
is the basic ID Hubbard model. On the right-hand side of Eq. 



(3) 



(4) 



is the kinetic-energy operator, t is the first-neighbor transfer integral, and U is the on-site repulsion associated with 
electronic double occupation. The operator 



n j, T n j, I > 



(5) 



counts the number of electron doubly-occupied sites and the operator 

a = CT o- , (6) 

counts the number of spin cr electrons at site j. The operators g and Cj. CT which appear in the above equations are 
the spin a electron creation and annihilation operators at site j respectively. Moreover, on the right-hand side of Eq. 
\i c = 2/i, [i s — 2/ip-ff, Mo is the Bohr magneton, and 



s: 



i 



[N a - N] 



1 



[JV T - 7VJ . 



(7) 
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are the diagonal generators of the SU(2) ry-spin S c and spin S s algebras [isLliil respectively. On the right-hand side 
of Eqs. (J2J) and Q the number of lattice sites N a is even and large, N a /2 is odd, N — J2cr N a , and N a = J2j ™j,cr- 

There are N-^ spin-up electrons and N± spin-down electrons in the chain of N a sites and with lattice constant a 
associated with the model Q). We assume periodic boundary conditions and employ units such that a = h = 1. When 



N a is odd the Fermi momenta are given by k F(J = ± 



where 



k Fa = 



7TN a 

N a 



or by kt a = 



(8) 



and 



When N a is even the Fermi momenta are given by kp a = kp a and k F(J = — kp a — jf- 

kp a = —kp a . Often we can ignore the 1/N a corrections of these expressions and consider k F<J ~ ±kFa = ±7rn CT and 
kp = [kp^ + kpi\/2 = Trn/2, where n a = N a /N a and n = N/N a . The electronic density reads n — rif + n±, and the 
spin density is given by m = — This paper uses the same notations as Ref. |3l| and often uses and refers to the 
equations introduced in that paper. For instance M. c ^±\li (and Af a ±1/2) denotes the number of ±1/2 holons (and 
±1/2 spinons), whereas L c> ±\ii (and £ s> ±i/2) denotes the number of ±1/2 Yang holons (and ±1/2 HL spinons). 

The Hamiltonian HsoU) defined in Eq. commutes with the six generators of the 77-spin S c and spin S s algebras 
|4H l49j , their off-diagonal generators reading, 



.1)3 J J 



3 



(9) 



and 



The Bethe-ansatz solvability of the ID Hubbard model is restricted to the Hilbcrt subspace spanned by the lowest- 
weight states (LWSs) of the r\ spin and spin algebras, i.e. such that S a — — |50| . 
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III. CPHS-ENSEMBLE SPACES AND GROUND-STATE DISTRIBUTIONS AND RAPIDITIES 

The studies of Ref. |3l| considered Hilbert subspaces associated with fixed values of 77-spin S c , spin S s , holon 
number M c , and spinon number M s . That was the choice suitable for the investigations on 77-spin-holon, spin-spinon, 
and charge c pseudoparticle separation of the excitations of the ID Hubbard model at all energy scales. However, 
fixing S c , S s , M c , M s does not fix the iVf and N\ electron numbers. 

In later sections of this paper and in Refs. |32l l37| transitions from a ground state with given spin a electron 
numbers to excited states are studied. Thus it is convenient for these studies to consider subspaces with fixed number 
of spin a electrons. From the use of Eq. JJJ and Eq. (35) of Ref. [211] > the N\ and N\ electron numbers can be 
expressed in terms of the M c +1/2 and M c -1/2 holon numbers and M s< _|_i/ 2 and M s ^ -1/2 spinon numbers as follows, 

A^T = 1{ N «~ ( M c,+i/2 - M c ,_ 1/2 ) + (M a , +1/a - M ai _ V2 )} , (11) 

and 

N.i = \{N a - (M Ci+1/2 - M c ,_ 1/2 ) - (M a , +1/2 - M a ,_ 1/2 )} , (12) 

respectively. We find from the use of Eq. (32) of Ref. |3l| that the M C)+1 / 2 holon and M s ,+1/2 spinon numbers are 
exclusive functions of the N c c-pseudoparticle number and M C) _i/ 2 holon and M s _ 1 / 2 spinon numbers and can be 
written as, 

Af Ci+1/2 = N a - N c - A/ c ,_ 1/2 ; M St+1/2 = N c -M 8i _ 1/2 . (13) 

It is convenient for our studies to express the and N± electron numbers in terms of the N c c pseudoparticle 
numbers, Mc,— 1/2 holon numbers, and M Si _i/ 2 spinon numbers, with the result, 

N^=N C + M Ci _ 1/2 - M Sj _ 1/2 , (14) 

and 

Ar; = M C) _ 1/2 + M a ,_ l/2) (15) 

respectively. 

We call an electron ensemble space a Hilbert subspace spanned by all states with fixed values for the and 
electron numbers. Furthermore, we call a CPHS ensemble space, where CPHS stands for c pseudoparticle, —1/2 
holon, and —1/2 spinon, a Hilbert subspace spanned by all states with fixed values for the numbers iV c , M c ^ _i/ 2 , and 
Mg t -i/2 of c pseudoparticles, —1/2 holons, and —1/2 spinons respectively. It follows from Eqs. and (fT5f) that 
in general, an electron ensemble space contains several CPHS ensemble spaces. This means that different choices of 
the numbers N c , M c _ 1 / 2 , and M S) _ 1 / 2 can lead to the same electron numbers and N^. In contrast, according to 
Eqs. I|14|) and l|15fl . a given choice of N c , M Cj _i/ 2 , and M S) _i/ 2 always corresponds to an unique choice for the Nf 
and N{ values of the electron numbers. 

Each CPHS ensemble space contains several subspaces with different numbers L c _ 1 / 2 of —1/2 Yang holons, L s _ 1 / 2 
of — 1/2 HL spinons, and sets of numbers {N c>v } and {N Si v } of c, v and s, v pseudoparticles respectively, corresponding 
to different i/ = 1,2, 3, ... branches. According to the results of Ref. |3J, the fixed values of each subspace obey the 
constraints imposed by the following equations, 



M Cj _ 1/2 = L Cj _ 1/ 2 + y £ d vN c , v ; M s ,_ 1/2 = L S) _ 1/2 + ^i/jV S;1 ,. (16) 

We call CPHS ensemble subspace a Hilbert subspace spanned by all states with fixed values for the numbers N c , 
L c _i/2, and L s _i/ 2 and for the sets of numbers {N c u } and {N Si „} corresponding to v — 1,2,3,... branches. It 
follows from Eqs. (|16f) that in general, a CPHS ensemble space contains several CPHS ensemble subspaces. This means 
that different choices of the values for L c _jy 2 and L s< _i/ 2 and sets of numbers {iV CjJ/ } and {N Si v } corresponding to 
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u = 1,2, 3, ... branches can lead to the same values for M Ci _i/2 and M s -1/2- hi contrast, a given choice of the values 
of N c , L c _ 1 / 2 , and L s _ 1 / 2 and sets {-/V C: „} and {N s ^ v } corresponding to v = 1,2,3, ... branches always corresponds 
to an unique choice of the values of N c , M C) -1/2, and M S; -1/2- 

Finally, we note that one does not need to provide the values of A/ c , +1/2 and M s +1 / 2 in order to specify a CPHS 
ensemble space. These numbers are given by Eq. (|13f) and are not independent. Furthermore, from combination of 
Eqs. (30), (32), and (34) of Ref. [3l| one finds that 

OO 

L c ,+i/2 = N a - N C -2J2" N c>v - L Ci _i /2 ; 

!/=l 

OO 

L.,+1/2 = N c -2^vN s , v -L s> _ l/2 . (17) 

!/=l 

Thus L c + i/2 and £ 8) +i/2 are not independent and one does not need to provide these values in order to specify a 
CPHS ensemble subspace. Therefore, often we do not consider the values of the holon numbers M C)+ i/2 and L c , +1/2 
and of the spinon numbers M s 4-1/9 and £ S) 4-1/2 in the expressions considered below. 

According to the studies of Ref. [3l|, our choice of using the numbers N c and M c _i/ 2 to label the charge sector of a 
CPHS ensemble space corresponds to a description of charge transport in terms of electrons. In this case the elementary 
charge carriers are the N c chargeons and M c 1/2 —1/2 holons, which carry charge — e and — 2e respectively. Here 
— e denotes the charge of the electron. The alternative use of the numbers iV c and M c+ i/ 2 to define the charge sector 
of a CPHS ensemble space would correspond to a description of charge transport in terms of electronic holes. In this 
case the elementary charge carriers would be the N c rotated electronic holes and M c +1 / 2 +1/2 holons, which carry 
charge +e and +2e respectively j-jll . Also the transport of spin has two alternative descriptions which correspond to 
the —1/2 spinons and +1/2 spinons respectively. Here we have chosen to label the spin sector of the CPHS ensemble 
spaces by the number M Sj _ 1 / 2 of —1/2 spinons. 

The finite-energy Landau liquid functional discussed in Ref. [3l| and used in the ensuing sections corresponds to the 
CPHS ensemble subspace description in terms of the values of the numbers of —1/2 Yang holons and —1/2 HL spinons 
and of the values of the numbers of c pseudoparticles at band-momentum q and of a, v pseudoparticles belonging 
to a = c, s and v = 1,2, 3, ... branches at band-momentum q. According to the studies of Ref. [31|, such a choice 
corresponds again to the description of charge transport in terms of electrons. In this case the c, v pseudoparticle 
carries charge — 2ve. 

IV. ELECTRON DOUBLE OCCUPATION OF THE PSEUDOPARTICLES, YANG HOLONS, AND HL 

SPINONS 

Let us denote the number of doubly occupied sites by D. By number of doubly occupied sites we mean here the 
expectation value of the operator (0) , 

D={D)=J2(4 A c j , t cl l c jA ). (18) 
3 

Another useful quantity is the kinetic energy T, which we define as the expectation value of the operator flljl. 

T = ^ XM a Cj+1, c) + (4 + 1, a Cj, a}} ■ (19) 

3: ° 

Let |^) be an energy eigenstate of the Hamiltonian © of energy Eh — (ip\Hjj\il>) given by the general expression, 

N 

E H = -2t-± / dqN c {q) cos k(q) 
2vr J q - 

N °° fit, v ( I -i 

+ 4 ^E/_ dqN e , v {q)Re{^l-{k c , v {q)+ivU/4£)^ + UL c> _ 1/2 , (20) 



6 



where k(q) and A c t/ (q) are the rapidity functionals defined by the integral Eqs. (13)-(16) of Ref. jU and N c (q) 
and N Cjll (q) are pseudoparticle band- momentum distribution functions also defined therein. The expressions of the 
band-momentum limiting values of the integrals on the right-hand side of Eq. 112011 are given in Eqs. (B14)-(B17) 
of Ref. l31l . This energy expression is obtained by combining the energy expression provided by the Bethe-ansatz 
solution |22ll38Ll4t| with the SO{A) symmetry of the model and the holon and spinon description introduced in Ref. 
|3lj . Note that such a general energy expression refers to the whole Hilbert space of the ID Hubbard model. The term 
involving the number of —1/2 Yang holons vanishes in the case of the LWSs of the SU(2) 77-spin algebra associated 
with the Bcthe-ansatz solution and is not provide d by such a solution. 

According to the Hellmann-Feynman theorem [3J,|35j, the electron double-occupation D, and kinetic energy T, 
associated with the above energy eigenstate are given by 



= W ^u w = au = ~mT ' (21) 



and 



rs|# | w=1 « sl |, (22) 

respectively, where Hjj is the Hamiltonian defined in Eq. |(3J and the energy Eh is given by Eq. H2U|) . 

The use of expression <|21|) in the evaluation of the electron double-occupation spectrum associated with creation 
of a c pseudoparticle, a, v pseudoparticle, —1/2 Yang holon, or —1/2 HL spinon involves the pseudoparticle Landau 
liquid energy functional. This functional was introduced in Ref. |3Sl ] and expressed in terms of holon and spinon 
numbers in Ref. l3ll. 



A. THE PSEUDOPARTICLE, HOLON, AND SPINON LANDAU-LIQUID ENERGY FUNCTIONAL 



The energy Landau-liquid functional introduced in Ref. |38j is generated by expanding the c and a, v pseudoparticle 
band-momentum distribution functions N c (q) and {N a ^(q)} where a = c, s and v = 1,2,3, .... Such an expansion 
is performed around the ground-state values of these functions given in Appendix C of Ref. |3lT | and inserting these 
distributions in the energy functional E associated with the Hamiltonian ijTjl. Such an energy functional is given by 



E = E H + J2 ^S« + - [L D> +1/2 - L c> _ 1/2 - -±] , (23) 

a— c, s 

where En is the functional l|2(J|l . the parameters \x a are the same as in Eq. |QJ, and S™ are the eigenvalues of the 
diagonal generators defined in Eq. (7J|. Let us consider that the initial ground state corresponds to values of the 
electronic density such that < n < 1 and values of the spin density such that < m < n. Throughout this paper 
we consider such ground states. The studies of Ref. ^38] have not expressed the energy Landau-liquid functional in 
terms of the deviations in the value of the —1/2 Yang holon and —1/2 HL spinon numbers. That problem was first 
addressed in Ref. [3ll | . As discussed in that reference, the general excitation spectrum of interest for the problem of 
the finite-energy correlation functions can be written as the sum of a finite energy loq and of a gapless contribution 
expressed in terms of the pseudoparticle energy bands. Such a energy Landau- liquid functional reads, 



AE = w + -M dqe c (q)AN c (q) + -^ / dq e a , 1 (q) AN S , 1 (q) 

J Qr — Q.t. 1 



'9c J-qs,i 



+ E E / dqel„(q)AN Q ^(q) 7 (24) 

a=c, s u=l+5 a s 1<*- " 



where AN c (q), AN c , L ,(q) 1 and AN sl/ (q) are the pseudoparticle band-momentum distribution function deviations. 
The band-momentum limiting values q s ,i, q a ,v, and q^r are defined in Eqs. (B14)-(B17) of Ref. |3l| . whereas the 
pseudoparticle energy bands e c (q), e Si i(q), and e° v (q) are defined by Eqs. (C15)-(C21) of the same reference. 

The expression for the energy parameter ujo on the right-hand side of Eq. (|24|l which controls the finite-energy 
physics is determined by the number of —1/2 holons and number of —1/2 spinons (except those associated with s, 1 
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pseudoparticles) of the final states relative to the ground state. Creation of a —1/2 holon requires an amount of 
energy 2/i and momentum 7r whereas creation of a —1/2 spinon (except those which are part of s, 1 pseudoparticles) 
requires an amount of energy 2/io-ff and zero momentum. Thus the energy parameter loq is a simple expression in 
terms of the deviations in the numbers of —1/2 holons and —1/2 spinons which reads, 



uj = 2/i AM, _ 1/2 + 2p H [AM,, _i/ 2 - AN S , J , (25) 

where AM Q , _ 1 / 2 = [^ Q ,-i/2 — M° _ x / 2 ] are the deviations in the numbers of —1/2 holons (a — c) and of —1/2 
spinons (a = s) and AN Si i = [N s , i — N® x ] is the deviation from the ground state value of the number of s, 1 
pseudoparticles. The ground-state pseudoparticle numbers are given in Appendix C of Ref. [3l| . 

The final states associated with the energy functional 1)24(1 belong to different CPHS ensemble spaces, depending 
on the numbers of quantum objects of these states. According to the form of the general energy spectrum defined in 
Eqs. (|24|l and (|25|) . final states belonging to CPHS ensemble spaces with finite occupancies of —1/2 holons have finite 
excitation energy relative to the ground state. These states belong to CPHS ensemble subspaces which are spanned 
by states with finite occupancies of —1/2 Yang holons and/or of c, v pseudoparticles. Again following expression l|25|l . 
for finite values of the spin density m, states belonging to CPHS ensemble subspaces spanned by states with finite 
occupancies of —1/2 HL spinons and/or of s, v pseudoparticles belonging to v > 1 branches, have also finite excitation 
energy relative to the ground state. However, most of the results of the present paper consider the limit of zero spin 
density, m = 0, where the energy spectrum of these states is gapless. 

On the other hand, for general values of n and m the low-energy Hilbert subspace corresponds to the CPHS 
ensemble space which contains the ground state. The corresponding ground-state c pseudoparticle, —1/2 holon, and 
— 1/2 spinon numbers are given in Eqs. (C24) and (C25) of Ref. Such a CPHS ensemble space is spanned by 

states with finite pseudoparticle occupancy for the c and s, 1 pseudoparticle bands only. We recall that according 
to the results of Ref. [31j the s, 1 pseudoparticle is a composite quantum object. It is constituted of two spinons of 
opposite spin projection. In the studies of Refs. |5lL l52l l53j the c and s, 1 pseudoparticle branches were referred to 
as c pseudoparticles and s pseudoparticles respectively. The excited states contained in that CPHS-cnscmble space 
have band-momentum distribution functions whose occupancies are different from the ones of the c and s, 1 ground- 
state distributions provided in Eqs. (CI) and (C2) respectively, of Ref. [IJ. At fixed spin a electron numbers, the 
occupancies of these excited states can be generated from the ground state by pseudoparticle - pseudoparticle hole 
processes in the c and s, 1 bands. 



B. THE GROUND-STATE ELECTRON DOUBLE OCCUPATION 

In the particular case of a zero spin-density ground state, from the use of the general Eq. I|21|l we find 

where Do denotes the ground-state electron double occupation and Eqsh stands for the ground-state energy, 

N f q Fc N f Q 

E GS h = -It— \ dq cosk°(q) = -2t— / dk 2np c (k) cos k . (27) 
2?r J q - c 2ir J_ Q 

Here the Fermi band-momenta qp c are defined by Eqs. (C4)-(C7) of Ref. j^H and the ground-state functions k°(q) 
and 2irp c (k) are studied below. This energy expression corresponds to the Hamiltonian 10 and is obtained by use 
of the ground-state numbers and deviations of Eqs. (C1)-(C3) and (C24)-(C25) of Ref. [31| in the general energy 
functional J5DJ). In order to study the ground-state double occupation (|26|l . let us provide useful information on the 
related energy (J37J. The ground-state function k°(q) on the right-hand side of Eq. I)27|) is related to the distribution 
2ixp c (k) considered below by the following equation, 



dk°{q) _ 1 
d( l ~ 27rpo (*"(?)) 



(28) 
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The function k°(q) is also related to the ground-state rapidity functions A° v (q) such that, 

<9A° Jq) 1 

T = 7 r; a = c,s, i/ = 1,2,3, .... (29) 

99 27rcr a)I/ (A0 iV (g)) 

The above functions 2np c (k) 1 27T(T c . v ( A) , and 27rcr S)1 /(A) are the solutions of the following integral equations, 

4icosfc U 



2 ^ = 1 + — ]_ B dK WTWU^W 2 ™ s ' l{A) ' (30) 

f At \ At f Q vU 

2na c „(A) = 2Re[ ■ =) / dk -—- ... , , —2np c {k), (31) 

, K) V( 4 *) 2 - {AtA-wUfl 7T J_ Q (^[/) 2 + (4i) 2 (sinfc-A) 2 Hcyj y ' 



and 



At f Q vU 

2ira s ^{K) = — dk , — 2np c (k) 

7r J_q (yU) 2 + (4i) 2 (smfc - A) 2 

- (^ i )2^,.(A'), (32) 

where the expression of the function Qj^sr) and the parameters Q and £> are given in Eqs. (C22) and (C23) of Ref. 
[3l| . We note that 2np c (k) and 27rcr Sl i(A) equal the functions 2izp(k) and 27rcr(A) respectively, introduced by Lieb 
and Wu j2l|. At half filling and spin density m — the parameters Q and B read Q — tt and B = oo and Eqs. 1(30(1 
- (1321) can be solved by Fourier transform. Use of the obtained expressions in Eq. 1(28(1 and in Eq. 1(29(1 for a, v = s, 1 
leads to the following expressions, 

poo sinl x sin k°(q) ) 

q = k (q) + 2l dxJ (x) ^ — y ; [ ? |<7r, (33) 

and 



o x [1 + e 2t 



/■oo sin(xA° x (g) 

g = / dx J (x) V ' / ; \q\ < it 12 , (34) 
Jo x cosh(^) 

where Jo(x) and Ji(x) are Bessel functions. These expressions define the inverse functions of k°(q) and A" 
respectively. Based on these expressions it is straightforward to find that at half filling and zero spin density the 
energy l|2"7Jl can be written in closed form with the result |21| . 

E GSH = -AN a tr dx J J^^. (35) 
Jo x[l + e 2t J 

One can use the obtained ground-state energy expressions in the evaluation of the electron ground-state double 
occupation 1(26(1 . We find that in the limits of vanishing and infinite on-site repulsion U/t the function f(n, U/t) of 
Eq. 1)26(1 is given by 1 and respectively. Closed form expressions for that function can be obtained for electronic 
densities < n < 1 in the limits U/t — > and U/t >> 1, with the result, 

f(n,U/t) = 1; U/t->0, 

/At\ 2 / sin (27m J. 
f(n,U/t) = (-) ln2(l- ^ ; ); U/t»l. (36) 
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For the specific case of half filling, n = 1, one finds from the use of Eqs. (|2(j[) and (|35[l the following closed- form 
expression for the function /(n, U/t) 35], 



(i 



J (x) Ji(x) 
l + cosh(^) 



According to the general expression (|2*2|) . the ground-state kinetic energy T is given by T = tdEcsn/dt [34]], 
where Egsh stands again for the ground-state energy associated with the Hamiltonian 

C. THE ELECTRON DOUBLE OCCUPATION FUNCTIONAL 

Our goal is the study of the deviations AD from the ground-state electron double occupation Do given in Eq. I|26|l . 
generated by changes in the numbers of pseudoparticles, —1/2 Yang holons, and —1/2 HL spinons associated with 
transitions to other energy eigenstates. Such a study provides information on the U/t dependence of the electron 
site distribution configurations of the real-space lattice which describe the energy eigenstates. The above-mentioned 
deviations are interesting quantities which provide useful information about the localization/delocalization degree of 
the corresponding elementary quantum objects. Fortunately, the electron double-occupation D of the final energy 
eigenstates can be obtained by means of the general expression (|21|l . Thus one can compute the electron double 
occupation deviations AD = [D — Do]. These can be obtained simply by taking the U derivative of the excitation 
energy defined by the following Landau-liquid energy functional, 

AE H = AE- A*aA5?-y[AL Ci+1/a -AL c ,_ 1/2 ], (38) 

a— c, s 

where AE is the energy functional 1)24(1 ■ AL c ±i/ 2 denotes the deviations in the numbers of ±1/2 Yang holons, and 
the remaining quantities are the same as on the right-hand side of Eq. (|20|l . 

The deviations AD can be written in functional form. To first order in the —1/2 Yang holon and —1/2 HL 
spinon number deviations and pseudoparticle band-momentum distribution function deviations, the double-occupation 
functional reads, 

AD = J2 AL Q ,-i/ 2 D Q ,_ 1/2 + ^ T C dqAN c {q)Dc{q) 

a=c,s •'la 



9c . 



TV 

+ E| dqAN a , v {q)D a ,u{q) . (39) 

Here AL a _ x / 2 is the deviation from the ground-state number of —1/2 Yang holons (a — c) and —1/2 HL spinons 
(a = s) and AN c (q) and AN a .y(q) are the c and a, v pseudoparticle band-momentum distribution function deviations 
given in Eq. (58) of Ref. |2l| respectively. On the right-hand side of Eq. (|39|l D a ^x/2 is the —1/2 Yang holon 
(a = c) and —1/2 HL spinon (a = s) double occupation, while D c (q) and D atV (q) denote the c pseudoparticle 
and a, v pseudoparticle double occupation spectra respectively. These double occupations and spectra equal the 
corresponding deviation in the value of double occupation l|18(> which results from creation of a —1/2 Yang holon or 
— 1/2 HL spinon and of creation of a c or a, v pseudoparticle at band- momentum q respectively. The pseudoparticle 
double-occupation spectra can be expressed as, 

D c (q)= gu ; D a>v (q) = — ( 40 ) 

where the pseudoparticle bands e c (q) and e aj „(g) are given by, 

e c (<?) = e c (q) + U/2- fj, + fi H ; e s , i(q) = e s> i(q) - 2^ H ; 

tc.viq) = £ v (q)+vU; e Sl u = e° SiV (q); v>\. (41) 
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Here the pseudoparticle energy bands e c (q), e s , 1 (q), e° v (q), and v {q) for v > 1 are denned by Eqs. (C15)-(C21) 
ofRef. ml. These are the energy bands of the energy functional 124|) . The —1/2 Yang holon and —1/2 HL spinon 
electron double-occupation deviation numbers are simply given by 

£>c-i/2 = i; ^,-1/2 = 0. (42) 

As for electron double occupation, we could also introduce a kinetic-energy functional associated with the kinetic- 
energy deviations AT = [T — T ]. However, the only information we need for our present study is whether creation 
of an elementary quantum object leads to a finite or a vanishing kinetic-energy deviation AT. 

The pseudoparticle electron double occupation spectra l|40|) and —1/2 Yang holon and —1/2 HL spinon electron 
double occupation (|42|l provide interesting physical information about the degree of localization/delocalization of the 
corresponding quantum objects. For instance, Eq. I)42|) reveals that for the whole parameter space creation of a 

— 1/2 Yang holon leads to AD — 1. That excitation involves creation of two electrons. From the expression of the 
77-spin flip generators given in Eq. 0, it follows that creation of a —1/2 Yang holon involves the creation of an extra 
on-site electron pair on a ground-state empty site for all real-space lattice electron site distribution configurations 
which describe the initial state. Consistently, we find that creation of a —1/2 Yang holon leads to a kinetic-energy 
deviation AT = 0. Since creation of such a quantum objects leads to no change in the value of the ground-state 
kinetic energy, we say that the —1/2 Yang holon carries no kinetic energy. Therefore, in spite of having charge — 2e, 

— 1/2 Yang holons do not contribute directly to charge transport, as mentioned in Ref. [3l]]. This concept of an 
elementary quantum object carrying kinetic energy when its creation leads to a finite deviation AT, is inspired in the 
electronic conductivity sum rule. That sum rule refers to the frequency dependent electronic conductivity and states 
that it is proportional to |T| [23 . Thus if creation of a quantum object does not change the value of T, it follows 
that the conductivity sum rule remains unchanged. Therefore, we say that such a quantum object carries no kinetic 
energy. Furthermore, we find from analysis of Eq. I|42|) that creation of a —1/2 HL spinon leads to no change in the 
ground-state double occupation. This is consistent with the form of the generators defined in Eq. i|10[l . Thus the 
spin flip associated with creation of a —1/2 HL spinon corresponds to electrons located at singly occupied sites for all 
real-space lattice electron site distribution configurations which describe the initial ground state. 

Let us start the discussion of the pseudoparticle electron double-occupation spectra defined in Eq. (|4*U1> by consid- 
ering the limit U/t — > 00. These spectra can be easily evaluated and in that limit the elementary quantum objects 
have simpler expressions in terms of real-space lattice electron site distribution configurations. As discussed in Refs. 
[3TI |3^|. such U/t — > 00 configurations are the same which describe the energy eigenstates for all values of U/t in 
terms of an effective electronic lattice. Such a lattice refers to the occupancy configurations of the rotated electrons 
introduced in Ref. [22 • It has the same number of sites N a and lattice constant a as the real-space lattice [2H [12 • In 
the limit of U/t — > 00 electron double occupation D is a good quantum number which equals the number of rotated- 
electron doubly-occupied sites 0, [j^] • We recall that for rotated electrons double occupation is a good quantum 
number for all values of U/t. As a result we find that as U/t — > 00 the electron double-occupation functional 
simply reads, 

AT 00 fla, v 

AD = AM c ,_ 1/2 = AL c ,_ 1/2 + -^J2 dqvAN CtV {q); U/t -► 00 . (43) 

From comparision of Eqs. (|39[1 and (|43|l . we arrive to the following U/t — ► 00 expressions for the pseudoparticle 
electron double-occupation spectra, 

D c (q)=D. tV (q)=0; D c . v {q) = u- U/t -> 00 . (44) 

These electron double-occupation spectra expressions show that increasing the number of c pseudoparticles by one does 
not change electron double occupation when U/t — > 00. On the other hand, from Eq. (|44|) we find that increasing 
the number of c, v pseudoparticles by one leads to an increase v in electron double occupation D. Moreover, we 
also studied the corresponding deviations in kinetic energy and found that in this limit the c pseudoparticles (c, v 
pseudoparticles) carry kinetic energy (carry no kinetic energy). 

These results reveal that in the present limit, the c pseudoparticles correspond to electron singly occupied sites and 
that the —1/2 holons and +1/2 holons correspond to electron pairs and electronic hole pairs at the same real-space 
lattice site respectively. In reference [3^ it is found that such a property is closely related to the rotated-electron 
effective electronic lattice mentioned above. As we discuss below, we find that at finite values of U/t a composite 
2i/-holon c, v pseudoparticle with band-momentum q such that \q\ < [it — 2k p] carries kinetic energy. This means 
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that for finite values of U/t the 2v extra electrons are not always created at the same site and thus that composite 
quantum object contributes to charge transport. However, we find that AT is a decreasing function of U/t which 
vanishes in the limit U/t — > oo. Thus creation of a c, v pseudoparticle leads in that limit to AT = and these objects 
become localized and carry no kinetic energy. It follows that in spite of having charge —2 ue, for U/t — > oo the c, v 
pseudoparticles do not contribute directly to charge transport, as mentioned in Ref. 31]. In that limit they describe 
a number v of localized electron pairs, each pair localized at the same site. 

On the other hand, in agreement with the results of Ref. |3lj | the chargeons have charge — e and we find that in the 
U/t — » oo limit they carry kinetic energy and thus are the only carriers of charge. When U/t — » oo the finite-energy 
Hilbert subspace is spanned by states with D — because it requires an infinite energy, D U, to excite the system 
into states with D > electron double occupation. In this limit, the c pseudoparticles become the spinless fermions 
of Refs. 155L loll IS?} • The associated chargeons describe the charge excitations associated with the electron singly 
occupied sites of the real-space lattice. The spin degrees of freedom of the electron singly occupied sites are described 
by the spinons. In the U/t — ► 00 limit creation of ±1/2 spinons describes spin flips of localized electrons at singly 
occupied sites. Therefore, creation of spinons does not contribute to double occupation, as confirmed by the results 
of Eq. Q43JI for the 2i/-spinon composite s, v pseudoparticles. This holds true for all finite values of U/t in the case of 
the HL spinons, as confirmed by Eq. 1 (4*21 . 

D. ELECTRON DOUBLE OCCUPATION SPECTRA FOR FINITE VALUES OF U/t AND ZERO SPIN 

DENSITY 

Let us now study the pseudoparticle electron double-occupation spectra of Eq. (|40|l for finite values of U/t. For 
simplicity, let us again consider the case of zero-spin density m — 0. In that case we find the following relations, which 
are valid for values of electronic density in the domain < n < 1 and for all values of U ft, 



The electron double-occupation spectra are derived from the expressions of Eq. I|41l) for the energy bands e c (q) 
and e c , „(<?) by means of the relation provided in Eq. (|40|) . Therefore, the features of these bands fully control and 
determine the band-momentum q, electronic density n, and on-site repulsion U ft dependence of the pseudoparticle 
electron double-occupation spectra. As we discuss in the ensuing section, the energy width of the band e c (q) is 
independent of U/t, being constant for all values of U/t and n and is given by W c = |e c (7r) — e c (0)| = At. In 
contrast, the width W c ^ u = |e c ,„(7r — 2kp) — e Cj „(0)| of the c, v pseudoparticle energy bands is a decreasing function 
of U/t, vanishing as U/t — > 00. The same behavior occurs also for the s, 1 pseudoparticle band, which is such that 
W s , 1 = \e s ,!(k F ) - e Si i(0) I -> as U/t ^00. 

This qualitative difference between the c pseudoparticle and the remaining pseudoparticle energy bands leads 
to qualitative differences in the corresponding electron double-occupation spectra as well. For instance, there are 
qualitative differences between the q dependence of the electron double-occupation spectrum D c (q) on the one hand, 
and that of the electron double-occupation spectra D St i(q) and D cu (q) on the other hand. While the function D c (q) 
changes slowly with band- momentum q, for finite values of U/t the functions D s ^i(q) and D CtV {q) show a significant 
q dependence. For instance, the values D s \(kp) = and D Ctl/ (jr — 2kp) = v are U/t independent, whereas the 
values of D s i(0) and D Cj „(0) depend on this quantity. In contrast, the U/t dependence of the spectrum D c (q) 
is very similar for all band-momentum values q. For finite values of U/t the electron double-occupation spectrum 
Ds.i(q) = Ds,i(—q) is a decreasing function of \q\, having its minimum value D s> i(kp) — 0, at \q\ = q s ,i = kp and 
its maximum value at \q\ = 0. For finite values of U/t the electron double-occupation spectrum D CjV {q) — D CtV (—q) 
is an increasing function of \q\, having its minimum value at \q\ = and its maximum value D c u (tt — 2kp) — v, at 
M = q c ,y = [tt- 2k F \. 

Since both the electron double occupation spectra D St \(q) and D c u (q) are not dependent on U/t for the limiting 
band- momentum values q = ±kp and q = ±[ir — 2kp] respectively, and show the strongest U /t dependence for q = 0, 
below we consider mainly the U/t and n dependences of the q = parameters Z? Si i(0) and -D Cj „(0). In agreement 
with the inequalities of Eq. (|46|l . for increasing values of the band- momentum absolute value \q\, these spectra reach 
smoothly the [//^-independent limiting values D s .i(±kp) = and D c _ u (±\-k — 2kp}) = v respectively. According to 
the results obtained in Ref. [32|, at m — only the c pseudoparticle, s, 1 pseudoparticle, c, 1 pseudoparticle, and 



DsAl) =0, 



v > 1 , 



(45) 



D C , V (Q) < D c ^{q) <D c , u {±[k -2k F ]) = u; 
D Sil (±k F ) = < D,,i(?) <D,,i(0). 



(46) 
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c, 2 pseudoparticle bands contribute significantly to the energy spectra of one-electron and two-electron excitations. 
Thus, below we limit our studies of the spectrum £) Cj „(0) to the branches v = 1 and v = 2. In the case of the 
D c {q) spectrum, we consider the n and U/t dependence of both the parameters D c (0) and Z? c (— 2fcp) — D c (2kp). 
Although the spectrum D c (q) is gently dependent on q, in addition to the value at the band-momentum q = we also 
consider the value at ±2k F . On the one hand, we study the parameter D c (2kp) because creation and annihilation of 
c pseudoparticles at the Fermi points ±2kp are associated with important one-electron and two-electron excitations 
pTj . On the other hand, since for the other pseudoparticle branches we study the U/t and n dependence of the 
electron double-occupation spectrum at q = 0, for completeness we also consider in our study the parameter D c (0). 

In figures 1 and 2 the electron double-occupation pseudoparticle parameters D c (0) and D c (2kp) are plotted as 
a function of U jt for different values of the electronic density n. The electron double-occupation pseudoparticle 
parameters D Si i(0), D c i(0), and D C:2 (0) are also plotted as a function of U/t and for the same values of the 
electronic density in Figs. 3-5. The curves plotted in these figures show that D c (0) and D c (2kp) are decreasing 
functions of U/t. Both of these functions have the maximum value D c (0) = D c {2kp) — n/2, for U/t — > and vanish 
in the limit U/t — ► oo. On the other hand, D Si i(0) vanishes both for U/t — ► and U/t — > oo and has its n dependent 
maximum value for a value of U /t between 1 and 4. While for small values of the electronic density n the parameters 
Dc, i(0) and -D c , 2(0) are decreasing functions of U/t, for values of n close to 1 there is a minimum value at a small but 
finite value of U/t. These parameters are given by n and 2n respectively for U/t — ► 0, and tend to 1 and 2 respectively 
as U/t — > 00. 

Before discussing the physical information contained in the curves plotted in Figs. 1-5, we consider some limiting 
cases where we could evaluate analytical expressions for the pseudoparticle electron double-occupation spectra. In 
addition, before that discussion we also indicate the cases where creation of the elementary quantum objects leads to 
a vanishing kinetic-energy deviation AT. In the limit U/t — ► we find the following expressions valid for densities 
< n < 1, 



D D (fi) = D D (±2k P ) = n/2; £> M (0) = 0; D c , v (Q) = vn. 

On the other hand, in the limit U/t » 1 the same quantities can be written as follows, 



(47) 



D c (0) 
D c (±2k F ) 

D c , u (0) 



2t\2 
U 

-2t\ 2 



n ln(2) 



sin(27rn) 



2irn 



(zt \ r 
—J n ln(2) |^2sin 2 (7rn) + 1 - 



sin(27rn) 



2irn 



2t^ 2 7m [ sin(27rn) 



C// 2 
4<\2 1 



UJ 2v 



2-KTl 

(1-n)- 



sin(27r(l - n)) 



2tt 



(48) 



Moreover, the following expressions are valid for all values of U/t > in the limit of vanishing density n — * 0, 



£> c (0) = £> c (±2fc F ) =£> s ,i (0) = 0; 



In the case of density n = 1 we could derive the following expressions, which are valid for all values of U/t > 



(49) 



P OO 

D c (0) = L» c (±2fc F ) = / dx- 



Ji(x) 
+ cosh(^) ' 

1 f 00 sinhf — ) 

D.,i(0) = - / dxJr^)— ttJ4; £» c ,„(0)= l ,. (50) 

2 7 [cosh(^)] 2 

Comparison of the expressions H47 [) -I|5U [) with the curves of the Figs. 1-5, reveals there is agreement between these 
expressions and the curves. 

As we mentioned above, we can also evaluate the kinetic energy deviations AT = [T— Tq] associated with elementary 
excitations. In the present study we are mostly interested in whether the deviations AT which result from creation of 
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a pseudoparticle, —1/2 Yang holon, or —1/2 HL spinon are finite or vanish. Importantly, by combining our analysis of 
the AT deviation expression with the related properties of the electron double-occupation spectra, we find that when 
such a deviation vanishes, creation of these quantum objects is either associated with creation of localized electrons 
or involves on-site spin-flip processes. From our study of the kinetic energy (|19|) we find that the deviation AT which 
results from creation of pseudoparticles, —1/2 Yang holons, or —1/2 HL spinons vanishes in the following cases: 

(i) Creation of one or a vanishing densitiy of s, 1 pseudoparticles at band momentum values q = ±/cpj for all values 
of U /t, n, and m. 

(ii) Creation of one or a vanishing densitiy of c, v pseudoparticles at band momentum values q = ±[tt — 2kp] for all 
values of U/t, n, and m. 

(iii) Creation of one or a vanishing densitiy of s, v pseudoparticles belonging to branches such that v > 1 at band 
momentum values q = ±[fci?t — fejpi] for all values of U/t, n, and m. 

(iv) Creation of —1/2 Yang holons for all values of U/t, n, and m. 

(v) Creation of —1/2 HL spinons for all values of U/t, n, and m. 

(vi) Creation of a, v pseudoparticles with a = c, s and v — 1,2, 3, ... for m = 0, all band-momentum values q, all 
densities n, and U/t — > oo. 

E. ELECTRON DOUBLE OCCUPATION SPECTRA AND THE PSEUDOPARTICLE DEGREE OF 

LOCALIZATION /DELOC ALIZ ATION 

We close this section with the analysis and a discussion of the information provided by the dependence on the elec- 
tronic density n, on-site electronic repulsion U/t, and band-momentum q of the electron double occupation quantities 
defined in Eqs. 14(il) - (|50|l and Figs. 1-5. We find that the dependence of the electron double-occupation quantities and 
spectra on these parameters provides insight into the degree of localization/delocalization of the —1/2 Yang holons, 
— 1/2 HL spinons, and pseudoparticles. 

While only spin a electrons can be created into the many-electron system or annihilated from it, creation of spin a 
electrons generates back-flow effects which can change the numbers and the occupancy configurations of pseudopar- 
ticles, holons, and spinons |32| . In the absence of electron addition or removal these elementary quantum objects are 
confined within the many-electron system, their occupancy configurations describing the exact energy eigenstates. As 
a result of the non-perturbative character of the electronic correlations, the description of the pseudoparticles, holons, 
and spinons in terms of electron site distribution configurations of the real-space lattice is a very complex many-body 
problem |3l| . 

For finite values of U/t the electron double occupation D defined in Eq. I|18f) is not a good quantum number. This 
just means that the ground state and the other energy eigenstates can be expressed as a superposition of real-space 
lattice electron site distribution configurations with different numbers of electron doubly occupied sites. On the other 
hand, in the limit U/t — > oo both electron double occupation D, double unoccupation, spin-down single occupation, 
and spin- up single occupation become good quantum numbers (3^, y^] • Thus in this limit all energy eigenstates can 
be expressed as a superposition of real-space lattice electron site distribution configurations with the same value of 
electron double occupation D. For instance, in the particular case of the ground state and finite-energy excited states 
that value is D = 0. Thus, the ground state electron double occupation Dq is such that Dq — > as U/t — > oo, as 
confirmed by the form of the expression defined in Eqs. i|26|) and (|36|l . 

It is useful for our discussion to express the deviations AiVj and AN^ in the electron numbers in terms of the 
deviations in the c pseudoparticle, a, v pseudoparticle, —1/2 Yang holon, and —1/2 HL spinon numbers. Based on 
the form of the expressions given in Eqs. (|14(l - (|16|l . we find the following relations, 



oo oo 



ATVf = AiV c - AL S , _ 1/a + AL C , - l/2 + ^v AN C , v -^v AN S , 



S, V 1 



(51) 



V=\ I/=l 



and 



oo oo 



AN V = AL S , _ 1/2 + AL C , _ 1/2 + v an c, v + Y." ANs > 



(52) 



!/=l V = l 



We recall that the numbers of +1/2 Yang holons and of +1/2 HL spinons are dependent on the numbers of the 
remaining quantum objects. In particular, it follows from Eq. (|17f) that the deviations AT Cj+1 /2 and AL s + i/ 2 can 
be expressed as follows, 
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OO 

AL c . +1/2 = -AN c -2j2v&N c>v - AL Ci _ 1/2 ; 
i/=i 

oo 

AL Si+1/2 = AN c -2j2"&N St „-AL Si _ 1/2 . (53) 

For simplicity, let us assume again that the initial ground state has zero-spin density m — 0. In order to study 
the degree of localization/delocalization of the quantum objects, we consider elementary excitations associated with 
creation of c pseudoparticles, s, 1 pseudoparticles, c, v pseudoparticles, and —1/2 Yang holons. According to Eqs. 
(|24H and l|25|) . for initial ground states with finite spin density m in the range < m < n the energy spectrum of 
excitations involving creation of s, v pseudoparticles belonging to v > 1 branches and/or —1/2 Yang holons is gapped. 
The study of the electron double occupation of such finite-spin-density excitations is an interesting problem which 
will be considered elsewhere. 

We start by considering the creation of a c pseudoparticle at band-momentum q. Creation of a c pseudoparticle 
is associated with creation of a spin- up electron, as confirmed by Eqs. I|51|) and 152|) . From Eq. I|53[l we find that 
such an excitation also involves annihilation of a +1/2 Yang holon and creation of a +1/2 HL spinon. According to 
Eqs. 1(43(1 and 1)48(1 . in the limit U/t — > oo this elementary excitation leads to no change in the ground-state electron 
double occupation Dq, Eq. l|2tj|) . Therefore, in this limit creation of a c pseudoparticle corresponds to creation of 
an extra electron singly occupied site and annihilation of an empty side of the real-space lattice. For decreasing 
values of U/t, the U/t dependence of the c pseudoparticle electron double occupation of Figs. 1 and 2 reveals that 
this excitation leads to a small positive finite electron double-occupation deviation AD = D c (q). At both q = 
and q = ±2kp this deviation is maximum when U/t — > 0, where from Eq. (|47|l and Figs. 1 and 2 we find that 
AD = D c (0) = D c (±2kp) — n/2. The deviation AD = D c (q) changes slowly with band-momentum q for all values of 
U/t and electronic density n. We emphasize that the kinetic-energy deviation AT which results from this excitation 
is always finite. We thus conclude that the c pseudoparticle does not describe localized electron configurations. 
The c pseudoparticle has it highest derealization degree in the limit of U/t — > oo, where it leads to no electron 
double-occupation deviation. In this case it becomes the non-interacting spin- less fermion of Refs. [{3EHEHE3- 

Let us now consider the case of the s, 1 pseudoparticle. Creation of a s, 1 pseudoparticle at band-momentum q is 
associated with a spin flip which creates a spin-down electron and annihilates a spin-up electron. This is confirmed 
by analysis of Eqs. (|51|l and 1(52(1. Moreover, it follows from Eq. I(53|l that creation of a s, 1 pseudoparticle also 
involves annihilation of two +1/2 HL spinons. The corresponding electron double occupation deviation is given by 
AD = D Si i(q) and equals the s, 1 double-occupation spectrum. According to Eqs. ((46(1 . (|47|l . and l(48|) . this spectrum 
has its minimum value at the band- momentum values q = ±kp, where it vanishes. In addition, we find that at this 
band-momentum value the s, 1 pseudoparticle does not carry kinetic energy. Thus at q — ±fcp and for all finite 
values of U/t and n creation of the s, 1 pseudoparticle is associated with an on-site electronic spin flip process in the 
real-space lattice. It follows that at this band-momentum values the s, 1 pseudoparticle has a localized character. For 
decreasing values of \q\ and finite values of U/t, the value of D Sj i(g) increases and we find that the s, 1 pseudoparticle 
carries kinetic energy. It follows that for |g| < hp, the spin flip electronic process has no pure on-site character. 
The delocalization degree of the s, 1 pseudoparticle is maximum for band-momentum q — 0. In contrast to the c 
pseudoparticles and c, v pseudoparticles, at finite values of U ft the degree of delocalization of the s, 1 pseudoparticle 
is in general highest when the value of D s \(q) is maximum. This can be understood by considering the U/t — -> oo 
limit. In this case the q = ±kp physics is extended to all values of q and the s, 1 pseudoparticle carries no kinetic 
energy. Moreover, according to Eqs. 1(43(1 and 148fl the electron double-occupation spectrum is in this limit given 
by D Si i(q) — 0. As was mentioned before, creation of the s, 1 pseudoparticle at any value of q for U/t — > oo is 
associated with an on-site electronic spin-flip process. It follows that in the case of the s, 1 pseudoparticle the highest 
degree of localization is achieved when it describes electron singly-occupied site excitations only. Thus, as the value 
of the on-site repulsion U/t is decreased and AD = D s i(q) becomes finite, the degree of delocalization increases. 
This effect is stronger for low values of band-momentum q. For the limiting values q = ±kp the s, 1 pseudoparticle 
describes a localized excitation for all values of U/t and n, as mentioned above. The curves plotted in Fig. 3 show 
that AD = D Sj i (0) is maximum for an intermediate value oiU/t. Interestingly, in the limit of U/t — > one finds that 
D s , i(<?) = for all values of q, as in the limit U/t — > oo. However, in this case the s, 1 pseudoparticle carries finite 
kinetic energy for all band-momentum values q except at q = ±kp. This means that for vanishing values of U/t the 
delocalization degree is not directly related to electron double-occupation. 

Another important elementary quantum object is the —1/2 Yang holon. Creation of a —1/2 Yang holon is associated 
with creation of two electrons of opposite spin projection, as is confirmed by Eqs. (|51|l and H52fl . From Eq. (|53|> 
we find that creation of a —1/2 Yang holon implies annihilation of a +1/2 Yang holon. The —1/2 Yang holon does 
not carry kinetic energy and from Eq. I)42|) we find that its creation leads to a electron double-occupation deviation 
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value AD = 1. Therefore, creation of a —1/2 Yang holon is associated for all finite values of U/t with the creation of 
a new localized electron doubly occupied site and the annihilation of an empty site in all real-space lattice electron 
site distribution configurations of the initial ground state, as mentioned previously. This is fully confirmed by the 
expression in terms of electronic operators of the fy-spin flip generators given in Eq. Jj5J • 

Let us now consider the c, v pseudoparticle. This is a s c = composite quantum object which results from the 
combination of v —1/2 holons and v +1/2 holons [3l]]. Since each —1/2 holon is made out of two rotated electrons of 
opposite spin projection [8lH82l ]. it follows that creation of a c, v pseudoparticle at band-momentum q is associated with 
creation of a number v of spin-down electrons and creation of an equal number of spin-up electrons. (We note that the 
number of rotated electrons equals that of electrons.) This is confirmed by Eqs. (|51|l and (|52|l . According to Eq. 153f) . 
creation of a c, v pseudoparticle also implies annihilation of 2v +1/2 Yang holons. The electron double-occupation 
deviation originated by creation of a c, v pseudoparticle is given by AD — D c u (q). According to Eqs. i|46[l . (|47|l . and 
(I48|l and Figs. 4 and 5, this deviation has its maximum value for the band-momentum values q = ±[tt — 2kp], where 
it is given by v. In addition, at this band-momentum values the c, v pseudoparticle does not carry kinetic energy. It 
follows that at these values of q creation of the c, v pseudoparticle is associated with creation of v localized electron 
double occupied sites and annihilation of v electron empty sites in all real-space lattice electron site distribution 
configurations of the initial ground state. Note that this result is valid for all electronic densities n and finite values of 
the on-site repulsion U/t. For these band-momentum values the c, v pseudoparticle corresponds to localized real-space 
lattice electron site distribution configurations similar to the ones obtained by creation of v —1/2 Yang holons. For 
decreasing values of \q\ and finite values of U/t the value D c u (q) decreases and we find that the c, v pseudoparticle 
carries finite kinetic energy. Thus, in contrast to the case of the s, 1 pseudoparticle, in the present case the degree of 
derealization increases with the decreasing of the electron double-occupation deviation value. For finite values of U/t 
and band-momentum values \q\ < [n — 2kp] the 2is created electrons are characterized by some degree of derealization 
which is maximum at q = 0. This means that for these band-momentum values the 2v electrons have some degree 
of derealization and are not all created in pairs at the same site of the real-space lattice. On the other hand, as 
U/t — »■ oo the q = ±[tt — 2kp] physics is extended to all values of q. In this limit the c, v pseudoparticle does not carry 
kinetic energy for all values of q and according to Eq. (|48|l the electron double-occupation spectrum is in that limit 
given by D CtU (q) — v. As discussed above, in this limit creation of a c, v pseudoparticle for all values of q is associated 
with creation of v localized electron double occupied sites and annihilation of v electron empty sites in all real-space 
lattice electron site distribution configurations of the initial ground state. Finally we emphasize that according to 
the total-momentum expression (36) of Ref. 31], the c, v pseudoparticle momentum spectrum is (n — q). Thus the 
band-momentum q = is associated with a contribution to the momentum given by (it — q) = it. This reveals that 
the maximum degree of derealization of the c, v pseudoparticles corresponds to a finite excitation momentum of 7r. 
Note also that for all the remaining pseudoparticle branches the momentum spectrum equals the band-momentum q. 

F. U/t -> oo ELECTRON DOUBLE-OCCUPATION SELECTION RULES 

Let the operator Oj^ be a product of a finite number 

la, h=±l 

of one-electron creation and/or annihilation operators. Here J\f/N a is vanishing small in the present thermodynamic 
limit and Ni c ,i 3 is the number of electronic creation and annihilation operators for l c = — 1 and l c = +1 respectively, 
and with spin down and spin up for l s = — 1 and l s — +1 respectively. 

Let us consider states Oj\f\GS) which result from application of a general A/"-electron operator Otf onto a ground 
state \GS) or onto any eigenstate of the spin a electron numbers. The pseudoparticle, holon, and spinon numbers 
of the ground state are given in Eqs. (C24) and (C25) of Ref. 31]. We consider transformations generated by 
application of a AA-electron operator O^f on any eigenstate of the spin a electron numbers and not necessarily on the 
ground state. Let \4>) be such a general state. It is not required that such an arbitrary state be an energy eigenstate. 
The corresponding state Oj\f\4>) is also an eigenstate of the spin a =f , J, electron number operators. It belongs to an 
electron ensemble space whose electron numbers differ from the numbers of the initial state by deviations AiVj and 
AN i , such that AN = AN^ + A Ni . These deviations can be expressed in terms of the numbers J\fi c ,i e of Eq. (|5"4"|) as 
follows, 



AiV t = E H)4+i; AN i = E Hc)-M c ,-i • (55) 

la = ±l l c = ±l 
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Following the results of Ref. [3l|, let us consider the four expectation values R a ,l a ~ (Ra,i a ), where a = c,s and 
l a = — 1, +1 and the corresponding operator i? Cj _i counts the number of electron doubly-occupied sites, R c . +i counts 
the number of electron empty sites, i? s . -i counts the number of spin-down electron singly-occupied sites, and R s ,+i 
counts the number of spin-up electron singly-occupied sites. The expressions of these operators are given in Eqs. (23) 
and (24) of Ref. [3l| . These operators obey the following relations 

Rc,+i = N a -N + Rc,-i; 4,-1 = ^-4,-1; = #r - K-i. (56) 

From the summation of the three relations of Eq. I|56|) wc find the following fourth dependent relation 

Ki a =N a . (57) 

a—c.s l a =±l 

The operational relations l|56(l - (|57(l are valid for the whole parameter space. Furthermore, they reveal that at given 
spin a electron numbers, out of the four expectation values of the operators of Eq. I|57f) onr y one is independent. In 
the previous section we have studied the double-occupation expectation value D = i? c> _i. The use of the operational 
Eqs. I|56|l - (|57|l provides the corresponding values for R c + i, R s _i, and R s ,+i- 

In the limit U/t — > oo the energy-eigenstate expectation values of the four number operators R c . _i, R c . +i, R s _ _i, 
and R s ,+i become good quantum numbers [sil l32j . Thus in that limit each energy eigenstate corresponds to real- 
space lattice electron site distribution configurations with the same values for R c ,~i, R c ,+i, R s ,-i, and R s ,+i- For 
finite values of U/t each energy eigenstate corresponds in general to a superposition of real-space lattice electron site 
distribution configurations with different numbers of electron doubly-occupied sites, electron empty sites, spin-down 
electron singly-occupied sites, and spin-up electron singly-occupied sites. 

Let us use the relations of Eq. I|56|l to evaluate expressions for the expectation values R c . +i, R Sl _i, and i? s . + i for 
the particular case of a ground state, from the corresponding m = electron double-occupation expression defined by 
Eqs. I|26l) and (|36|) . We find that for m — and for all values of U/t and for values of the electronic density in the 
domain < n < 1, these ground-state expectation values read, 

R°c + i = N a - N + | g) f(n, U/t) , (58) 

and 

N 



-ft, _i — n e 



.,+1 " 2 



(^)f(n,U/t)], (59) 



respectively. Here f(n,U/t) is the function given in Eqs. H36J) an d JHZJl. Note that these ground-state expectation 
values change from jR° _ t = D = f (f Y R° c . +1 = [N a - N] + f (§) , and jR° _ t = R° s< +1 = f 1 - (f ) as U/t 

to _ x = Dq = 0, R® +1 = [N a — N], and R® _ 1 — i?° +1 = ^ when U/t — » oo and these quantities become good 
quantum numbers. 

Importantly, in the limit C//£ — > oo the four deviations Ai2 a , i a — [R a , i a i ] such that a = c, s and Z Q = — 1, +1 
of any eigenstate of the spin a =f, | electron number operators which result from application onto that state of the 
above general operator O^f are restricted to the following ranges, 



M+hl* < AD = &Rc,-i < M-Ll. , (60) 

!»=±1 l,=±l 
- £ ^ Ai? ^+1 ^ E > ( 61 ) 



Z„=±l f»=±l 



l c ,l e = ±l i c ,Z 3 = ±l 



(62) 
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and 

- Si c ,i m M u ,i.<AR ai+1 < J2 s i«-i.M c ,i., (63) 

i ,i B =±i i c ,i s =±i 

respectively. Obviously, given one of the four ranges of values denned by the inequalities (|60J) - H63J) . the other three 
follow from the operational relations of Eq. 1|56H • 

The maximum (and minimum) value of the inequalities and (|61|l (and I|tj2|) and (ESJ) are reached when all 
the electrons created by the operator Oj\f, in a number of = ±iAf—i,i„, transform an equal number of electron 
singly-occupied sites into electron doubly-occupied sites. In addition, it is required that all the annihilation operators 
of that operator, in a number of _ ±1 -A/+i, i s , transform an equal number of electron singly-occupied sites into 
electron empty sites for all the real-space lattice electron site distribution configurations of the initial state with a 
number of electron singly occupied sites larger than or equal to Af. Note that if the initial state includes distribution 
configurations whose number of electron singly occupied sites is smaller than J\f, then the maximum value of AD is 
smaller than J^i =±i-^-i.L an d thus the inequalities l (5Uj l and (|5T|l (and ijfi^ and JOHJ) are also satisfied. 

The minimum (and maximum) value of the inequalities (|6L)fl and <|tjl|) (and l|t)2|) and (ESJ) are reached when all the 
electrons created by the operator Oj\f, in number of , ]\[-\j s , transform an equal number of electron empty 

sites into electron singly-occupied sites. Furthermore, it is required that all the J^i =±i-^+i, h annihilation operators 
of that operator transform an equal number of electron doubly-occupied sites into electron singly-occupied sites for 
all the real-space lattice electron site distribution configurations of the initial state with a number of electron empty 
sites equal to or larger than =±1 N-i,i B and a number of electron doubly-occupied sites equal to or larger than 
Si =±i ^+1, la ■ Note that if the initial state includes distribution configurations whose number of electron empty sites 
is smaller than _ ±1 A/"-i, i g and/or whose number of electron doubly-occupied sites is smaller than =±1 A/+i, i 3 , 
the minimum value of AD is larger than — J^i =±i-^+i,2. ancl thus the inequalities ijSOjl and llHT|) [and l |62"jl and |JB3J] 
are also satisfied. 

It is found in Ref. that the U/t — > oo electron double-occupation selection rules associated with the inequalities 
(|60|l - Ht)3|l are related to the occurrence of other selection rules which concern the range of the deviations in the number 
of ±1/2 holons and ±1/2 spinons. The holon and spinon selection rules are valid for all values of U/t and play an 
important role in the evaluation of correlation functions |32l l36l 1371 ] . These selection rules are exact in the case of 
rotated-electron operators. In the case of the corresponding electron operators such rules provide the numbers of 
holons and spinons of the final states which contribute most significantly to the few-electron correlation functions 
|32|. This relationship between the U/t — * oo electron double-occupation deviations and the deviations in the number 
±1/2 holons and ±1/2 spinons follows from the unitary transformation which maps electrons onto rotated electrons. 
This transformation also maps the electron double-occupation operator onto an operator that counts the number of 
— 1/2 holons. The latter operator is nothing but the rotated-electron double occupation operator [32^. 

V. THE PSEUDOPARTICLE ENERGY BANDS 

In section IV we found that the pseudoparticle energy bands defined by Eqs. (C15)-(C21) of Ref. [U] play an 
important role in the electron double occupation spectra of the pseudoparticles, as confirmed by Eq. J3IJ. Moreover, 
in reference |37| it is found that the shape of the lines in the frequency/energy and momentum plane where the 
peaks or edges in the one-electron and two-electron spectral-weight distributions are located corresponds to the 
pseudoparticle energy bands. Thus the study of the band-momentum dependence of these energy bands is a problem 
of interest for the understanding of the one-electron and two-electron spectral properties of the quantum liquid. The 
c pseudoparticle and s, 1 pseudoparticle energy bands were studied and plotted in Ref. |5l| for finite values of the spin 
density. On the other hand, the c, v pseudoparticle energy bands and s, v pseudoparticle energy bands for v > 1 were 
first introduced in Ref. |38||. yet these bands were not plotted in that reference. Since most studies on the finite-energy 
spectral properties of low-dimensional materials refer to zero magnetization, it is worthwhile plotting and discussing 
the band-momentum dependence of the above pseudoparticle bands for zero spin density. 

The general excitation energy spectrum of interest for the problem of the finite-energy correlation functions is given 
in Eq. I|24l) . The band- momentum independent term ojq is defined in Eq. (|25|l . The value of this finite energy is 
determined by the holon and spinon deviations and correspond to the finite-energy edges of the correlation functions. 
The other term of the energy spectrum (|24|l is of gapless character and in the case of zero magnetization involves 
the pseudoparticle e nerg y bands e c (q), e St i(q), and e° v (q) where v — 1,2,3,.... These bands are defined in Eqs. 
(C15)-(C17) of Ref. \3$. As discussed below, at zero magnetization both the energy and band-momentum width of 
the energy bands e° v {q) vanishes in the case of the v > 1 branches. According to Eqs. (C1)-(C3) of Ref. jU, in 
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the case of the zero-spin-density ground states the band e c (q) is filled (and empty) for band- momentum values such 
that \q\ < 2kp (and 2kp < \q\ < 7r). On the other hand, for such ground states the band e St x(q) is full and the bands 
e c,u{l) where v — 1,2,3,... are empty. The expressions of the bands e c (q), e s .i(q), e° ^g), and e° 2 (g) are plotted 
in Figs. 6-9 as a function of the band-momentum q for spin-density m — and different values of U/t and of the 
electronic density n. 

As is confirmed by analysis of the energy-band curves plotted in Figs. 6 (a) and 6 (b), both the energy and band- 
momentum width of the c pseudoparticle bands are density, spin-density, and U/t independent and given by At and 
2ir respectively. From the U/t dependence of the curves of these figures, we find that the main effect of increasing the 
on-site repulsion U/t is the increasing of the energy-band width of the ground-state filled sea and the corresponding 
decreasing of the band-energy width of the corresponding unfilled region. The main effect of changing the electronic 
density is on the ground-state c pseudoparticle Fermi points, which are given by ±2kp. 

In contrast to the energy and band-momentum widths of the c pseudoparticles, the energy width of the a, v 
pseudoparticle bands plotted in Figs 7-9 is a decreasing function of U/t, vanishing in the limit U/t — > oo. This property 
is a result of the fact that all 77-spin and spin configurations are degenerated in that limit. The band- momentum width 
of the c, v pseudoparticle bands and of the s, v pseudoparticle bands belonging to branches such that v > 1 are density 
and spin-density dependent respectively, and given by 2[ir — 2/cf] and 2[kp] — kp^] respectively. Note also that the 
a, v pseudoparticle energy bands e° = e S) ±(q) — 2fiQ H and e° a y {q) vanish at the band-momentum limiting values, 
i.e. 

e° atV (±q a , v ) = 0, (64) 

where according to Eq. (C12)-(C14) of Ref. 31], g Ci „ = [ir — 2kp], g S] 1 = kp^, and q Sj v = [kp^ — kpi\ for the v > 1 
branches. 

For finite values of the spin density m and magnetic field H , the s, v pseudoparticle energy bands e° v (fj) have for 
the v > 1 branches both finite momentum width, given by the above expression 2[fcj?t — kp-/\, and finite energy with, 
given by |e° „(0)|. The energy band value at q = is negative, i.e. e° „(0) < 0. As the spin density m — > 0, both 
widths 2[kp^ — kp^] — ► and \e° s „(0)| — > vanish. Thus, for the v > 1 branches the s, v pseudoparticle energy bands 
collapse into the point {q = , e° ^(0) = 0} in that limit. Since Figs. 6-9 correspond to m = 0, these bands are not 
plotted in these figures. The finite-magnetization case will be considered elsewhere. 

The c, v pseudoparticle bands plotted in Figs. 8 (a)-(b) and 9 (a)-(b) have their maximum value at q — 0. These 
bands are thus inverted. As discussed below, as the electronic density approaches 1, both the energy and momentum 
widths of these bands vanish. Except for the case of the s, 1 band of Fig. 7, all remaining bands plotted in Figs. 
6-9 show for finite values of U/t a quadratic dependence on q in the vicinity of the band-momentum limiting values 
q = ±<7 C = ±7r and q = ±q c , v = ±[tt — 2kp\. In contrast, note that the s, 1 band has a linear q dependence in the 
vicinity of the band-momentum limiting values q = ±q s ,i = ikp. On the other hand, for finite values of U/t all 
pseudoparticle bands show a quadratic dependence on the band-momentum q in the vicinity of the point q = 0, as 
confirmed by the curves of Figs. 6-9. 

From analysis of Figs. 8 and 9, we find that the e° ^q) and e° 2 (q) bands vanish at the band- momentum limiting 
values q = iq c ,v = ±[tt — 2kp\. This behavior also occurs for finite values of the spin density m. For such a values 
also the e° v (q) bands vanish for q = ±q s , v = ±[fcF| — ^fj.]- The energy band e° v (q) has its maximum value at q = 0, 
with e° j,(0) > 0, as shown in Figs. 8 (a)-(b) and 9 (a)-(b). The energy band e° „(<?), for spin density m > has its 
minimum at q — 0, with e° ^(0) < 0, as in the case of the s, 1 energy band plotted in Fig. 7 for m = 0. 

In the Appendix we consider the limiting expressions of the pseudoparticle energy bands in the case of zero spin 
density, m = 0. Note that according to the expressions given in Eqs. ljAl|l - l|A3|l of the Appendix, the q dependence 
of the c pseudoparticle energy band in the vicinity of q = ±7r and of the c, v bands in the vicinity of q = is different 
for finite values of U/t and in the limit U/t — > 0. As mentioned above, for finite values of U/t the dependence on the 
band-momentum q is quadratic in the vicinity of these band-momentum values. Thus in that case the q derivative 
of these bands vanishes both at q = ±7r and q — respectively, as confirmed by the curves of Figs. 6-9 and by the 
expressions of Eqs. (|Alj) and (|A3(I of the Appendix. In contrast, analysis of the data of these figures and equations 
reveals that in the limit U/t —> the q derivative at these band-momentum values becomes finite and the q dependence 
in their vicinity is linear instead of quadratic. This singular change results from the non-perturbative character of 
the electronic correlations which leads to a different physics at U/t — and for U/t > 0. In the particular case of 
half- filling, this effect occurs for the c pseudoparticle band at the band-momentum Fermi points q = ±7r. In this case, 
that effect is associated with the singular Mott-Hubbard insulator - metal transition. For other electronic densities 
this effect occurs at unoccupied ground-state band-momentum values and therefore, is not associated with such a 
transition. The vanishing at half-filling of the band e° „(<?) of Eq. (|A6|I of the Appendix results from its collapse into 
the point {q = 0,e° „(0) = 0} as S — > 0, where S = (1 — n) is the doping concentration. This is confirmed by the 
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dependence on n (and 8) of expressions l|A3(l for that band. We emphasize that both the s, v bands for the v > 1 
branches and the c, v bands collapse into the point {q = , e° a v {§) = 0} with a = s and a — c respectively, as the 
spin density m — > and the doping 8 = (1 — n) — > respectively 

The terms of the Landau energy functional given in Eq. (I24|) describe the energy spectra of all one-electron and two- 
electron elementary excitations of the ID Hubbard model (31]. Within such a functional description, each elementary 
excitation is simply described by specific values of the pseudoparticle band- momentum distribution function deviations 
defined in Eq. (58) of Ref. |3jJ. These deviations describe: (a) pseudoparticle - pseudoparticle hole processes in the c 
and s, 1 bands which conserve the pseudoparticle numbers; (b) creation or annihilation of c and s, 1 pseudoparticles; 
(c) creation of c, v pseudoparticles and of s, v pseudoparticles with v > 1; and (d) creation of —1/2 Yang holons and 
of —1/2 LH spinons. 

All the energy spectra obtained previously in the literature, as for instance the ones studied in the Refs. |4ll l42l liM 
can be expressed in terms of the pseudoparticle energy bands plotted for m — in Figs. 6-9. The authors 
of these studies used for each specific elementary excitation different forms for the two Lieb and Wu equations pH 
or/and for the thermodynamic Takahashi's equations [22 provided in Ref. 31]. On the other hand, insertion in the 
energy Landau-liquid functional H24JI of band-momentum distribution function deviations suitable to these excitations 
leads to the same energy spectra. In addition, the energy functional (|24|l provides the spectrum for all values of the 
energy of all excitations which contribute to the one-electron and two-electron physics, as mentioned above. 



VI. CONCLUDING REMARKS 



In this paper we studied the electron double-occupation deviations which result from creation of c pseudoparticles, 
a, v pseudoparticles, —1/2 Yang holons, and —1/2 HL spinons, where a — c, s and u = 1,2, 3, .... All energy eigenstates 
can be described in terms of occupancy configurations of these quantum objects which are related to electrons and 
rotated electrons in Refs. 0, [32 . Here we introduced an electron double-occupation functional whose coefficients 
are pseudoparticle band-momentum dependent double-occupation spectra. Our investigations provided interesting 
information about the electron site distribution configurations of the real-space lattice which describe the energy 
eigenstates. The study of these spectra provided important information on the localization/deslocalization degree of 
the real-space lattice electron site distribution configurations which describe the elementary quantum objects. We 
found that for some of the pseudoparticle branches such a degree of localization/delocalization is strongly dependent 
on the value of the band momentum. We also considered U/t — > oo selection rules, which limit the ranges of the 
electron double-occupation deviations that result from excitations generated by AA-electron physical operators. These 
rules are the starting point for finding related selection rules, which are investigated elsewhere [32 . The latter rules 
are valid for all values of U ft and play an important role in the study of finite-energy correlation functions. 

This paper is a first step in the use of the general holon and spinon scenario in the study of the finite-energy spectral 
properties of the ID Hubbard model. In addition to the introduction and study of the electron double-occupation 
functional, in this paper we studied the dependence on the band-momentum q, on-site repulsion U/t, and electronic 
density n of the pseudoparticle energy bands associated with such a functional. As was mentioned in this paper, these 
bands play a central role in the spectral properties of the model [32j . 

The holon and spinon description of the ID Hubbard model introduced in Ref. 31] and further studied here is 
used in Ref. [33] in the finding of further useful microscopic information on the spectral properties of the model. By 
combining the results obtained here with that information, one can evaluate expressions for finite-energy one-electron 
and two-electron correlation functions [^J [3(| [33 . That is a problem of major importance for the study of the 
unusual finite-ene rgy spectral properties observed in real quasi-one-dimensional materials GUMSEllSIEGllliEl El 
[Til IT2I IT3I Il4l IT5L [la. IrH [LsL Il9|. Eel . Elsewhere it is confirmed that the holon and spinon des crip tion of the energy 
eigenstates introduced in Ref. j3j] and further developed here and in the related articles [U H3, Ha> H3 i s suitable 
for the successful study of these finite-energy spectral properties. 
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APPENDIX A: LIMITING EXPRESSIONS OF THE PSEUDOPARTICLE ENERGY BANDS 



In this Appendix we present the limiting expressions of the pseudoparticle energy bands in the case of zero spin 
density, m — 0. In the limits U/t — > and U/t >> 1 one can evaluate closed form expressions for the bands e c (q), 
e s, i(q), and e„ u (q) where a — c, s and v = 1, 2, 3, defined in Appendix C of Ref. |31j. The simplest case is that 
of the s,v energy bands for the v > 1 branches, which at zero spin density m — are given for all values of U/t 
and n by e° u (q) = for \q\ < (kp^ — kpi) =0. As m — > these v > 1 bands collapse in that limit into the point 
{q = , j,(0) = 0}. The expressions of the remaining pseudoparticle bands for zero spin density m = 0, values of 
the density < n < 1, and U/t — > and t//t >> 1 are the following, 



e c (g) 



-4t 
-2i 
-2f 

?7 



7rn\ 



cos (l) 

cosMgl 
cos(3) — cos(7rn) 



7rn\ /7ra\ 

t)- cos (t) 



g|<2/c F , £//*-► 0; 

2k F < \q\ < tt, U/t->0: 



2n ln(2) sin 2 ((j) - sin 2 (7m) 



kl < TT, cv* » 1 



(Al) 



1(9) 



-2t 

m 

U 



cos(g) 



7T 

2 



/7rn\ 

cos ItJ 

sin(27rn) 
2^ 



< fcf 



, <7 
cos l — 



U/t- 
\q\ < k F 



U/t » 1 



(A2) 



At cos 



17 



1 

2^ 



sin(27T(5) 



2tt 



<?| < (tt - 2fc F ) , U/t^0 

kl < (7T-2fcjr), f7/t» 1. 



cos 



(A3) 



where 5 = (1 — rz.) is the doping concentration. The expressions fA"Tj) and (£3) were already obtained in the c and 
s = s, 1 pseudoparticle studies of Ref. [HJ ■ From analysis of the band-momentum q, on-site repulsion U /t, and density 
n dependence of the band expressions ljAl(l - l|A3(l . we find that these agree with the corresponding dependence on the 
same parameters of the curves of Figs. 6-9. 

At half filling and spin density m = the integral equations which define the ground-state rapidities and pseu- 
doparticle bands can be solved by Fourier transform. Thus in this case one can derive the following expressions for 
the pseudoparticle bands, which are valid for all values of U/t > p2\ . 



e c (q) = -2t cos k°{q) + 1 



cos 

- At I dx Ji (x) 

'0 x \1 + 



(x sink°(q)j 



|9| <7T, 



(A4) 



e,,i(«) = -2t 



cos(xA° sl (q) 

dx J\ (x) 77— 

U xcosh(^) 



kl <7T/2, 



(A5) 



and 



4,u(q) = 0; \q\<(n-2k F ) = 0, (A6) 

The inverse functions of the half- filling ground-state rapidity- momentum function k°(q) and rapidity function A° 1 (q) 
involved in expressions (|A4ft and l|A5(l are given in Eqs. and (|3^| respectively. 
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FIGURE CAPTIONS 

Figure 1 - The electron double-occupation parameter D c (0) plotted as a function of U/t for different values of the 
electronic density n. 

Figure 2 - The electron double-occupation parameter D c (2kp) plotted as a function of U/t for different values of 
the electronic density n. 

Figure 3 - The electron double-occupation parameter D s ^ i (0) plotted as a function of U/t for different values of the 
electronic density n. 

Figure 4 - The electron double-occupation parameter D Ct i(0) plotted as a function of U/t for different values of the 
electronic density n. 

Figure 5 - The electron double-occupation parameter D Ct2 (0) plotted as a function of U/t for different values of the 
electronic density n. 

Figure 6 - The pseudoparticle energy band e c (q) in units of t plotted for electronic density (a) n = 1/2 and (b) 
n = 5/6 and on-site repulsion U/t -> 0, U/t = 2.0, U/t = 4.0, U/t = 10, U/t = 20, and U/t -> oo. The ground- 
state Fermi level corresponds to zero energy and is marked by a horizontal line which overlaps the e c (q) band at the 
band-momentum Fermi points q = ±2fcp = ±tt/2. 

Figure 7 - The pseudoparticle energy band e s . i(q) in units of t plotted for electronic densities n = 1/2 and n — 5/6 
and on-site repulsion U/t -> 0, U/t = 2.0, U/t = 4.0, U/t = 10, U/t = 20, and U/t -> oo. The ground-state Fermi 
level corresponds to zero energy overlaps the energy band at the band-momentum Fermi points q = ±kp. 

Figure 8 - The pseudoparticle energy bands (a) e° l (q) and (b) e° 2 (?) m units of t plotted for electronic density 
n = 1/2 and on-site repulsion U/t -» 0, U/t = 2.0, U/t = 4.0, U/t =10, and U/t = 20. The U/t -> oo limit of these 
bands corresponds to a horizontal line located at the zero energy level. 

Figure 9 - The pseudoparticle energy bands (a) e° 1 (g) and (b) e° c 2 (q) in units of t plotted for electronic density 
n = 5/6 and on-site repulsion U/t -> 0, U/t = 2.0, U/t = 4.0, U/t =10, and U/t = 20. The U/t -> oo limit of these 
bands corresponds to a horizontal line located at the zero energy level. 
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